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ABSTRACT 

In  wireless  networks,  relay  transmissions  can  enable  co¬ 
operative  diversity  by  forming  virtual  antenna  arrays.  The 
optimum  resource  allocation  in  such  relay  networks  is  crit¬ 
ical  to  enhance  their  performance  and  efficiency.  However, 
existing  works  on  resource  optimization  only  consider  single¬ 
relay  systems  and  focus  on  the  power  allocation.  In  this  pa¬ 
per,  we  consider  both  the  power  optimization  and  the  loca¬ 
tion  optimization  for  systems  with  arbitrary  number  of  relays. 
Equally  attractive  is  that  our  investigation  is  tailored  for  dif¬ 
ferential  modulations,  which  bypass  the  channel  estimation 
at  the  receiver  and  are  particularly  suitable  for  wireless  re¬ 
lay  networks.  We  first  derive  an  upper  bound  of  the  error 
performance.  Based  on  this  bound,  we  then  develop  the  op¬ 
timum  energy  and  distance  allocation  schemes  that  minimize 
the  average  system  error.  Analytical  and  simulated  compar¬ 
isons  confirm  that  the  optimized  systems  provide  consider¬ 
able  improvement  over  un-optimized  ones.  In  addition,  we 
show  that  location  optimization  may  be  more  critical  than  en¬ 
ergy  optimization. 

1.  INTRODUCTION 

Recent  studies  show  that  networks  consisting  of  unmanned 
robotic  and  tele-operated  aerial  and  ground  vehicles  etc.  serv¬ 
ing  as  sensors,  communications  relays  and  weapons  systems 
are  enablers  of  the  C4ISR  capabilities  [14].  In  particular,  re¬ 
lay  networks  provide  diversity  gains  by  forming  virtual  an¬ 
tenna  arrays  with  distributed  network  nodes  in  wireless  com¬ 
munications.  By  exploiting  the  cooperative  diversity,  the  an¬ 
tenna  packing  limitations  can  be  eliminated  and  the  spatial 
diversity  gain  can  be  achieved  [1,  7,  9].  A  majority  of  exist¬ 
ing  works  on  relay  networks  focuses  on  coherent  demodula¬ 
tion  based  on  the  availability  of  the  channel  state  information 
(CSI)  at  both  the  relays  and  the  destination  node  (see  e.g., 
[7,  9,  10]).  Accurate  estimation  of  the  CSI,  however,  can  in¬ 
duce  considerable  communication  overhead  and  transceiver 
complexity,  which  increase  with  the  number  of  relay  nodes 
employed.  In  addition,  CSI  estimation  may  not  be  feasible 
when  the  channel  is  rapidly  time-varying.  To  bypass  chan¬ 
nel  estimation,  cooperative  diversity  schemes  obviating  CSI 
have  been  recently  introduced.  These  relay  systems  rely  on 


noncoherent  or  differential  modulations,  including  conven¬ 
tional  frequency-shift  keying  (FSK)  and  differential  phase- 
shift  keying  (DPSK)  [3,  4,  13],  as  well  as  space-time  coding 
(STC)based  ones  [5,  11]. 

To  improve  the  error  performance  and  enhance  the  energy 
efficiency  of  relay  networks,  optimum  resource  allocation  re¬ 
cently  emerges  as  an  important  problem  attracting  increas¬ 
ing  research  interests  (see  e.g.,  [2,  8,  12]).  These  works  are 
based  on  different  relaying  protocols  (amplify-and-forward, 
decode-and-forward  and  block  Markov  coding),  under  var¬ 
ious  optimization  criteria  (signal-to-noise  ratio  (SNR)  gain, 
SNR  outage  probability,  energy  efficiency  and  capacity),  and 
with  different  levels  of  CSI  (instantaneous  CSI  and  channel 
statistics).  However,  all  of  them  only  consider  the  power  al¬ 
location  and  mostly  focus  on  a  single-relay  setup. 

In  this  paper,  we  consider  a  relay  network  with  arbitrary 
number  of  relays.  More  importantly,  we  treat  the  resource 
allocation  as  two  optimization  problems:  the  optimization  of 
the  energy  (power)  distribution  and  the  optimization  of  the 
relay  location.  Equally  attractive  is  that  our  analysis  is  tai¬ 
lored  for  relay  systems  with  differential  modulation,  which  is 
known  to  reduce  the  receiver  complexity  by  bypassing  chan¬ 
nel  estimation  [3,  4,  5,  13].  To  enable  the  resource  optimiza¬ 
tion,  we  first  derive  an  upper  bound  of  the  overall  symbol 
error  rate  (SER)  performance  for  relay  networks  employing 
the  decode-and-forward  (DF)  protocol.  The  energy  and  loca¬ 
tion  optimization  will  then  be  carried  out  based  on  this  perfor¬ 
mance  bound.  We  show  that  under  the  constraints  of  the  total 
energy  per  symbol  and  the  source  and  destination  distance, 
the  optimum  SER  performance  can  be  achieved  through  the 
joint  energy  and  location  optimization.  Interestingly,  location 
optimization  may  be  more  critical  than  energy  optimization. 
In  other  words,  near-optimum  performance  can  be  achieved 
by  location  optimization  alone,  but  not  by  the  energy  opti¬ 
mization  alone. 

The  rest  of  this  paper  is  organized  as  follows.  The  sys¬ 
tem  model,  including  the  relay  protocol,  the  differential  (de¬ 
modulation  and  the  diversity  combining  rules  for  the  relay 
transmission,  is  introduced  in  Section  2.  In  Section  3,  an 
SER  upper  bound  is  established  for  a  relay  setup  with  arbi- 
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source  s  relays  {rk}%=1 


Fig.  1.  Setup  of  the  relay  network. 

trary  number  of  relays.  The  optimum  energy  and  distance 
allocation  is  derived  in  Section  4.  Simulations,  comparisons 
and  discussions  are  presented  in  Section  5.  Summarizing  re¬ 
marks  are  given  in  Section  6. 

Notation:  We  use  (•)*  for  conjugate,  E[-]  for  expectation, 
■)?{•}  for  the  real  part,  and  :=  for  “is  defined  as”.  CM(y,  a2) 
represents  the  complex  Gaussian  distribution  with  mean  /i 
and  variance  a2. 

2.  SYSTEM  MODEL 


bypass  channel  estimation.  Specifically,  with  the  nth  phase- 
shift  keying  (PSK)  symbol  being  denoted  as  sn  =  ej27rc"/M, 
cn  G  {0, 1, M  —  1},  the  corresponding  transmitted  signal 
from  the  source  is: 


The  encoded  signal  is  broadcast  via  a  common  channel,  and 
the  received  signal  at  the  /cth  relay  is  given  by 

Vr„k’s  =  VW-X  +  zrn\  k  =  l,2,...,L,  (2) 

where  £s  is  the  energy  per  symbol  at  the  source  node,  the 
fading  coefficient  of  the  channel  between  s  and  rk  during  the 
nth  symbol  duration  is  ~  CM(0,  a2r  ),  and  the  noise 
component  z] [k  ~  CM(0,Mrk).  This  signal  is  differentially 
demodulated  and  remodulated  independently  at  each  relay  rk. 
The  demodulation  step  generates  an  estimate  from 
in  Eq.  (2),  using  the  decision  rule  that  we  will  present  in  the 
next  subsection.  The  remodulation  step  is  carried  out  as  in  Eq. 
(1),  but  with  sn  replaced  by  its  estimate  and  xsn  replaced  by 
.  Then,  the  received  signal  at  the  destination  correspond¬ 
ing  to  each  relay  node  is  given  by 

ydnrk  =  V^Kir*xrn*+zt  k  =  1,2, ..., L,  (3) 


Consider  a  network  setup  with  one  source  node  s,  L  re¬ 
lay  nodes  {rk}k-i  and  one  destination  node  d,  as  depicted 
in  Fig.  1.  Each  node  is  equipped  with  a  switch  that  controls 
its  transmit/receive  mode  to  enable  half-duplex  communica¬ 
tions.  Multiplexing  among  the  network  nodes  can  be  achieved 
via  frequency-division,  time-division  or  code-division  tech¬ 
niques.  For  notational  convenience,  we  will  consider  the  time- 
division  multiplexing  (TDM).  However,  the  presented  anal¬ 
ysis  and  results  are  readily  applicable  to  frequency-division 
multiplexing  (FDM)  and  code-division  multiplexing  (CDM). 

2.1.  Relaying  Protocol  and  Channel  Modeling 

We  consider  the  decode-and-forward  (DF)  relaying  pro¬ 
tocol,  in  which  the  relay  nodes  de-modulate  the  signal  from 
the  source  node,  then  re-modulate  and  forward  to  the  destina¬ 
tion  node.  During  the  first  segment  of  each  symbol  duration, 
the  source  node  broadcasts  the  first  symbol  to  all  relay  nodes. 
Next,  each  relay  transmits  the  remodulated  signal  to  the  des¬ 
tination  during  their  distinct  segments  within  the  rest  of  the 
symbol  duration.  As  a  result,  the  total  of  L  source-to-relay 
(s  —  r)  links  share  a  common  channel,  whereas  the  L  relay-to- 
destination  (r  —  d)  links  have  mutually  orthogonal  channels. 
With  TDM,  the  orthogonality  is  ensured  by  assigning  distinct 
time  segments  to  each  relay.  This  setup  is  widely  adopted  in 
relay  networks  (see  e.g.,  [2,  4,  9]). 

As  we  mentioned  in  the  introduction,  differential  modu¬ 
lation  is  considered  at  both  the  source  and  relay  nodes  to 


where  £Tk  is  the  energy  per  symbol  at  the  fcth  relay  node,  the 
fading  coefficient  of  the  channel  between  rk  and  d  during  the 
nth  symbol  duration  is  h^rk  ~  CM(Q,  crkd  ),  and  the  noise 
component  is  ~  CM (Q,Mrj).  Accordingly,  we  can  find  the 
received  instantaneous  signal-to-noise  ratio  (SNR)  between 
the  transmitter  j  and  the  receiver  i  as 

\h^\  %  ...  _ 

7 i,j  =  ^  ,  h3  e  {s,rk,d}. 

Then,  the  average  received  SNR  is  7 =  (<r2.  £j)/Mi. 

2.2.  Differential  Demodulation  and  Decision  Rules 

As  mentioned  before,  differential  demodulation  is  per¬ 
formed  at  the  relay  and  destination  nodes.  To  derive  the  de¬ 
modulation,  decision  and  diversity  combining  rules,  let  us  be¬ 
gin  with  the  received  signal  at  the  relay  or  the  destination  node 
Vn  =  hnxn  +  zn,  which  is  extracted  from  Eqs.  (2)  and  (3)  by 
dropping  the  superscripts.  Using  the  differential  encoding  in 
Eq.  (1),  the  received  signal  can  be  re-expressed  as: 

Vn  —  hn(.Xn—lSn)  “f“  Zn  =  yn—\Sn  ~\~  Zn  , 

where  z'n  =  zn  —  zn-\sn.  For  M-ary  PSK  symbols,  it  follows 
that  E[s*  sra]  =  1.  Hence,  the  conditional  distribution  of  yn  is 
complex  Gaussian  with  mean  yn-  i  sn  and  variance  2 A).  As 
a  result,  we  obtain  the  log  likelihood  function  (LLF)  of  yn  as: 

ImiVn)  :=  lnPy„|S„(yn|4n)  =  Vn-llm}  ,  (4) 


where  i,j£{s,  r k,  d},Im=e:’2nrn/M  and  m£{ 0, 1, M  —  1}. 

At  the  A:th  relay  node,  the  differential  demodulator  is  then 
straightforward: 

srnk=ej2nm'/M:m'  =  aTgmax  lrffs(yrnk's) 

m 

=  arg max  ^{{Vn^Tyn-Jm}  ■ 

m 

At  the  destination  node,  however,  there  are  L  different  LLF’s 
corresponding  to  the  L  transmitted  signals  from  the  relays: 

tirk(Vn)  =  M{(ydnrk)*  yt\Im}  ,  k  =  1,2, . . . ,  L  .  (5) 


If  the  channel  state  information  (full  or  partial)  is  known 
at  the  relays  and  the  destination  node,  then  it  is  possible  to 
combine  the  LLF’s  by  capturing  the  detection  error  at  the  re¬ 
lay  node  according  to  the  so-termed  transition  probability  (see 
e.g.,  [3]).  However,  keeping  in  mind  that  differential  modula¬ 
tion  is  considered  in  the  first  place  because  of  its  capability  of 
bypassing  channel  estimation,  we  will  focus  on  the  scenario 
where  no  channel  state  information  is  available.  Under  this 
circumstance,  the  LLF’s  in  Eq.  (5)  have  to  be  combined  with 
equal  weights.  Accordingly,  the  decision  rule  at  the  destina¬ 
tion  node  can  be  readily  obtained  as: 

L 

Sd=e*W/M  .  m'  =  argmax  J2^{(ydnrkrydn-\Im}  ■ 

k=  1 

With  no  channel  information  assumed  at  either  the  relays  or 
the  destination  node,  this  decision  rule  turns  out  to  be  the  dif¬ 
ferential  detection  with  postdetection  equal  gain  combining 
(EGC)  [17,  Chapter  6.6], 


3.  ERROR  PERFORMANCE  ANALYSIS 

To  facilitate  our  resource  optimization,  we  will  derive  the 
analytical  expression  of  the  error  performance  for  the  cooper¬ 
ative  system  described  in  the  preceding  section.  Symbol  error 
probability  of  cooperative  networks  with  relay  transmissions 
has  been  derived  in  [9]  for  coherent  detection,  and  in  [13] 
for  a  differential  scheme  with  a  single  relay,  both  employing 
the  amplify-and-forward  (AF)  protocol.  Here,  we  consider  a 
general  /.-relay  setup  under  the  DF  protocol,  with  differential 
(de-)modulation  and  diversity  combining  assuming  no  CSI. 


Let  us  denote  the  average  symbol  error  rate  (SER)  at  the 
fcth  relay  node  as  Pe,rk  ■  For  differential  M-ary  PSK  (DMPSK) 
signaling,  the  s  —  link  SER  Pe,rk  can  be  obtained  as  [16, 
Chapter  8.2.5] 


_ \/ 9 psk  l‘7T/2 M1b(-[1-  ^/l- gpsK  COS  0]) 

2n  J_n/ 2  1 —  \/l  —gpsK  cos  6 


(6) 


where  M7(x)  =  1/(1  —  27),  Vx  >  0,  and  7  represents  the 
average  signal-to-noise  ratio  (SNR).  In  particular,  for  M  =  2 
(DBPSK),  Eq.  (6)  can  be  simplified  as 


P, 


1 

2(1  +  7Vfc,s ) 


(7) 


At  the  destination,  the  signals  from  the  L  relays  are  com¬ 
bined  to  make  a  decision.  Conditioned  on  that  the  symbol  sn 
is  correctly  demodulated  and  remodulated  at  all  relay  nodes, 
the  conditional  SER  PBtd  can  be  obtained  by  applying  the  re¬ 
sults  for  /.-diversity  branch  reception  of  M-phase  signals  in 
[15,  Appendix  C]  as: 


Pe,d  = 


(— 1)L"1(1  -p2)L  (dL 


=cof 


<%L_1  \b-ff2  M 
! —pcos(: t/M) 


7t(L— 1)! 

psin(n/Af) 

\/b—p2cos2(n/M)  —  p2  cos2  (p/M)  J  J  /{,=] 


,(8) 


where  p  =  'yd,rk/(l  +  ld,rk )■  For  DBPSK,  Eq.  (8)  can  be 
simplified  as 


Using  the  unconditional  SER  Pe,rk  at  the  relays  and  the 
conditional  SER  PBid  at  the  destination,  we  formulate  an  up¬ 
per  bound  on  the  overall  average  error  performance,  namely 
the  unconditional  SER  Pe  at  the  destination,  as  follows: 


Proposition  1  With  Pe,rk  and  Pe  d  given  by  Eqs.  (6)  and  (8), 
respectively,  an  upper  bound  on  Pe  can  be  found  as: 

L 

Pe<Pe  =  l-l[(l~Pe,rk)(l-Pe,d).  (10) 

fc=  1 

Proof:  To  prove  that  Eq.  (10)  provides  an  upper  bound  on 
the  exact  SER  Pe,  let  us  start  with  the  probability  of  correct 
detection  Pc  =  1  —  Pe.  Counting  the  events  that  lead  to  the 
correct  detection,  Pc  can  be  obtained  as 

F(FPr{[K‘  =Sn,  Wfc)  n  (Sn  =  Sn)}  (li) 

u  [{srn^sn,  for  some  fc)  n  (s^  =  sn)]} 

=  Pr{ s^=sn | srnk=sn,\/rk}- Pr{ srf  =sn,  Vrfe } 

+  Pr{s((=sn|s(/^sn,for  some  fc}-Pr{s(/^sn,for  some  k}. 

where  Srf  and are  the  symbol  estimates  formed  at  the  relay 
/■/.  and  the  destination  d,  respectively.  The  first  summand  in 
Eq.  (11)  turns  out  to  be  IIfc=i(l  ~  Pe,Th){  1  -  Pe,d),  which 
leads  to  the  upper  bound  in  Eq.  (10). 

Several  remarks  are  due  here  on  the  second  summand  in 
Eq.  (11),  which  corresponds  to  the  gap  between  the  true  SER 
and  its  upper  bound  A Pe  :=  Pf.  -  If  and  determines  the 


Fig.  2.  SER  at  different  yT<a  values  (L  =  2,  M  =  2). 

tightness  of  the  error  bound  in  Proposition  1.  For  DBPSK 
with  a  single  relay  (L  =  1,  M  =  2),  this  gap  can  be  easily 
obtained  as: 

A  Pe  =  Pe,rPe,d ■  (12) 

For  practical  Pe  r  and  Pe^  values  (e.g.,  <  10-3),  APe  is  neg¬ 
ligible  compared  with  Pe  =  Pe  r  +  Ped  —  2 PerPe>d-  How¬ 
ever,  for  L  >  2,  all  possible  errors  have  to  be  considered  for 
both  the  s  —  r  and  r  —  d  links,  which  renders  APe  analyti¬ 
cally  untractable.  But  intuitively,  as  the  L  increases,  A Pe  also 
increases  since  there  is  an  increasing  chance  that  detection  er¬ 
rors  at  the  relay  nodes  do  not  lead  to  a  detection  error  at  the 
destination  node.  In  addition  to  this  effect,  the  performance 
bound  Pe  and  the  gap  APe  also  depends  on  the  quality  of  the 
s  —  r  and  r  —  d  links. 

These  effects  are  evident  from  the  simulated  examples  in 
Figs.  2  and  3,  where  a  relay  network  with  L  =  2  relay  nodes 
using  DBPSK  signaling  is  considered  at  various  7r  s  and  7 d,r 
levels.  In  these  simulated  examples,  the  channels  between  the 
source  and  all  relays  have  identical  powers  aft  =  <7^  , 

Vfc,  which  implies  that  yrk,s  =  7r,s>  Vfc.  Accordingly,  we 
have  Pe,rk  =  Pe,r>  Vfc,  and  Pe  <  Pe  =  1  —  (1  —  Pe^r)L{  1  — 
Pe,d)  from  Proposition  1.  Likewise,  the  SNR  between  all  the 
relay  nodes  and  the  destination  have  the  same  power,  Jd,rk  = 
Jd,r,  Vfc.  In  both  Figs.  2  and  3,  the  bound  Pe  closely  captures 
the  dependency  of  the  system  SER  on  the  SNR  levels  tv>s  and 
7 d,r-  Specifically,  we  have  the  following  observations: 

•  Fig.  2  reveals  that,  at  any  given  value  of  yr.s,  the  sys¬ 
tem  SER  exhibits  an  error  floor  as  7 d,r  increases.  Intu¬ 
itively,  this  error  floor  comes  from  the  detection  error  at 
the  relays,  which  heavily  relies  on  the  s  —  r  link  quality 
TV-^  and  can  only  be  reduced  by  imposing  sufficiently 
high  7r>s. 

•  On  the  contrary.  Fig.  3  shows  that,  at  medium-to-high 
7 d,r  levels,  the  overall  SER  can  always  be  reduced  by 


Fig.  3.  SER  at  different  7 d,r  values  (L  =  2,  M  =  2). 


increasing  the  SNR  of  the  s  —  r  link  7,%s  and  does  not 
exhibit  any  error  floor. 


4.  OPTIMUM  RESOURCE  ALLOCATION 

In  the  preceding  section,  we  derived  a  simple  upper  bound 
on  the  overall  SER  of  the  relay  network.  Our  simulated  and 
numerical  examples  indicate  that  the  s  —  r  link  SNR  and  the 
r  —  d  link  SNR  play  unbalanced  roles  in  determining  the  sys¬ 
tem  SER.  In  this  section,  we  will  investigate  the  effects  of 
resource  allocation  on  the  SER  performance.  We  will  show 
that  an  optimum  allocation  of  the  limited  resource  is  possible, 
and  it  achieves  the  optimum  system  error  performance. 

For  analytical  tractability,  we  consider  an  idealized  /.-relay 
system  with  all  relay  nodes  located  at  the  same  distance  from 
the  source  and  the  destination  nodes;  that  is,  Ds  ric  =  DSjr 
and  l)ri.  y  =  Dr.d,  Vfc.  It  is  then  reasonable  to  assign  equal 
energies  at  all  relay  nodes  £rk  =  £r,  Vfc.  To  carry  out  the 
optimization  in  the  ensuing  subsections,  we  will  also  make 
use  of  the  relationship  between  the  average  power  of  channel 
fading  coefficient  er3  and  the  inter-node  distance  l)hl  as 
follows: 

Oh^  =  c  '  Di hJ  e  {a,  r,  d}  ,  (13) 

where  v  is  the  path  loss  exponent  of  the  wireless  channel  and 
C  is  a  constant  which  we  henceforth  set  to  1  without  loss  of 
generality. 

4.1.  Optimum  energy  allocation 

Problem  Statement  1  For  any  given  source,  relay  and  des¬ 
tination  node  locations  (Ds,r  and  Dr  d,  or  equivalently  a ^ 
and  <j'fld  ),  and  the  total  energy  per  symbol  £,  determine  the 


optimum  energy  allocation  £s  and  £r  which  minimize  Pe  in 
Eq.  (10)  while  satisfying: 

L 

£S  +  J2^rk  =£s+L£r  =£  ■  (14) 

fc= 1 

Without  loss  of  generality,  assuming  that  all  noise  compo¬ 
nents  are  independent  and  identically  distributed  (i.i.d)  with 
J\frk  =  Md  =  A fo-  Then,  by  defining  the  total  SNR,  p  := 
£ /Wo,  the  transmit  SNR  at  the  source  node  ps  :=  £s/J\fo  and 
the  transmit  SNR  at  the  relay  nodes  pr  :=  £r/Afo,  the  energy 
constraint  can  be  re-expressed  as  the  SNR  constraint : 

p  =  ps  +  Lpr.  (15) 

Using  Eq.  (13),  the  average  received  SNR’s  at  the  relay  and 
destination  nodes  can  be  expressed  in  terms  of  the  transmit 
SNR’s  as: 


Fig.  4.  Optimum  energy  distribution  (p  =  lOdB,  L  =  1,  DBPSK). 


7r,s  =  Ps<?hra  =  PsDs,r 
and  7 d,r  =  pro\dr  =  prDfd  .  (16) 

As  a  result,  the  total  energy  constraint  (15)  can  be  further 
rewritten  as 


can  find  ps  as  in  Eq.  (18).  The  detailed  proof  of  the  exact  so¬ 
lution  and  its  simplification  for  the  optimum  energy  allocation 
can  be  found  in  [6].  From  Eq.  (18),  it  readily  follows  that  the 
energy  allocation  ratio  between  the  source  and  the  relay  nodes 
is 


P  =  lr,s/olr  s  +  Lyd,r/a2hd  r  =  7 r,sDusr  +  L^rDvr  d  .  (17) 

To  gain  some  insights,  we  start  from  a  single-relay  setup 
and  establish  the  following  result: 

Proposition  2  For  a  single-relay  setup  with  L  =  1,  at  given 
s  —  r  and  r  —  d  distances  Ds^r  and  Dr and  under  the  total 
energy  constraint  in  Eq.  (14),  the  optimum  energy  allocation 
£s  should  satisfy : 


£s  _  £DStry/2 

£r  \Dr,d) 


(23) 


Eq.  (23)  reveals  explicitly  that  the  optimum  energy  allocation 
heavily  hinges  upon  the  inter-node  distances.  In  addition,  the 
path  loss  exponent  of  the  wireless  channel,  v,  also  affects  the 
optimum  energy  allocation.  Interestingly,  the  £s/£r  ratio  is 
linear  in  Ds,r/ Dryd  only  when  v  =  2. 
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(18) 


and  correspondingly,  £ r  =  £  —  £s. 


This  solution  is  achieve  by  solving  the  first  order  condi¬ 
tions  under  medium-to-high  SNR  values.  Treating  the  SER 
bound  Pe  as  a  function  of  7r,s  and  yd,r,  we  have  the  first  or¬ 
der  conditions  for  the  optimum  solution 


Fig.  4  depicts  the  transmit  SNR  ps  obtained  from  the  op¬ 
timum  energy  allocation.  A  one-dimensional  setup  is  consid¬ 
ered;  that  is,  Ds>r  +  Dr  d  =  Ds  d-  The  system  parameters 
are:  p  =  lOdB,  L  =  1,  Ds  d  =  1  and  v  =  (1,2, 3, 4).  In 
Fig.  4,  the  simulated  optimum  ps  is  plotted,  together  with  the 
exact  analytical  value  in  Eq.  (22)  and  the  approximation  in 
Eq.  (18).  Although  the  approximate  expressions  Eq.  (18)  is 
obtained  under  the  high  SNR  assumption,  they  remain  very 
accurate  even  at  the  medium  SNR  of  lOdB. 


dPe 

dPe 

dfd,r 


-  A D%  =  0 

-  A DvrA  =  0 


P  -  (7 +  7 d,rD"  d)  =  0  , 


(19) 

(20) 
(21) 


where  A  is  the  Lagrange  multiplier.  From  Eqs.  (19)  and  (20), 
we  can  achieve  the  expression  of  fd ,r-  Then,  by  substituting 
7 d,r  into  Eq.  (21),  the  following  exact  solution  (see  Eq.  (22)) 
for  ps  is  achieved.  By  applying  high  SNR  approximation,  we 


From  Fig.  4,  we  also  observe  that,  for  all  v  values,  the 
source  node  energy  £s  increases  as  the  relay  moves  towards 
the  destination  node.  With  v  =  2,  £s  increases  linearly  with 
Ds  r.  At  higher  values  of  the  path  loss  exponent,  v  >  2,  we 
observe  that 

Ps  f  <  D 's,r/ Ps,dt  when  Ds  r  <Z  Ds  d /2  ,  (241 

p  \  >  DStr/Ds,d,  when  Ds>r  >  Ds,d/ 2  • 

In  other  words,  the  optimum  energy  allocation  favors  the  link 
with  larger  inter-node  distance.  When  the  path  loss  exponent 


(22) 
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i/  =  1,  Fig.  4  shows  the  opposite  of  Eq.  (24). 

So  far,  we  have  been  focusing  on  the  single-relay  case, 
where  an  analytical  solution  in  Eq.  (22)  can  be  obtained  and 
a  very  accurate  and  insightful  approximation  is  available  un¬ 
der  the  high  SNR  assumption.  For  L  >  2,  however,  the  first 
order  conditions  obtained  by  differentiating  the  SER  bound 
Pe  have  complicated  forms,  which  render  analytical  solutions 
impossible.  Fortunately,  the  SER  bound  Pe  as  in  Proposition 
1  still  allows  for  a  numerical  search,  as  opposed  to  Monte 
Carlo  simulations  needed  otherwise. 


For  example,  with  DBPSK  and  L  =  2  ,  we  have 
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and,  accordingly,  the  SER  bound  is  given  by 


(25) 


Fig.  5.  Optimum  location  distribution  (p  =  lOdB,  L  =  1,  DBPSK). 


5  _  1  7r,s (2  +  7d,r)(l  +  27d,r)2  ^ 

6  4(l  +  7r,S)2(l  +  7d,r)3  '  1  ^ 

By  using  the  first  order  conditions  in  Eq.  (19)  and  the  high 
SNR  approximation,  the  optimum  fT_s  and  rid.r  should  satisfy 


4(1  +  7d,r)(2  +  7d,r)(l  +  ^ld,r)  _  ^ r,d 
37r,s(l  +  7r,s)  Ds,r 

Although  an  analytical  solution  is  not  readily  available,  one 
can  resort  to  the  numerical  search. 


4.2.  Optimum  distance  allocation 

Consider  now  a  one-dimensional  relay  setup  with  the  dis¬ 
tance  between  the  source  and  the  destination  being  Dsd .  If 
the  transmit  energies  at  the  source  and  relays  are  preset,  where 
is  the  optimum  location  to  place  the  relays?  To  answer  this 
question,  we  treat  the  distance  Dsd  as  a  fixed  resource  and 
formulate  an  optimization  problem  as  follows: 

Problem  Statement  2  For  any  given  transmit  energies  at  the 
source  and  relay  nodes  (£s  and  £r,  or  equivalently  ps  and  pr), 
and  the  path  loss  exponent  v  of  the  wireless  channel,  deter¬ 
mine  the  optimal  location  of  the  relays,  Dsr,  which  minimizes 
Pe  in  Eq.  (10)  while  satisfying  0  <  DS}7.  <  DS}d  ■ 

Starting  with  the  single-relay  ( L  =  1)  setup  and  apply¬ 
ing  the  high-SNR  approximation,  we  establish  the  following 
result: 


Proposition  3  For  a  single-relay  setup  with  L  =  1  and  the 
source-destination  distance  Ds  d,  and  let  £s  and  £r  denote 
the  prescribed  transmit  energy  levels  at  the  source  and  relay 
nodes,  respectively,  the  optimum  location  of  the  relay  is 

i/(— i) 

Ds,r  =  1/(*  1/(y_1}  •  Ds,d  ,  (28) 

Ps  +  Pr 

and,  accordingly,  Dr  d  =  Ds  d  —  DS}r. 


Similar  to  the  energy  optimization  case,  the  solution  can 
be  achieved  by  solving  the  first  order  conditions,  in  which 
we  treat  the  SER  bound  Pe  as  a  function  of  Ds  r  and  Dr.d- 
We  omit  detailed  proof  due  to  the  space  limit.  With  the  one¬ 
dimensional  setup.  Proposition  3  can  be  represented  as 


Dr,d  \Pr)  \£r ) 


(29) 


Interestingly,  Eq.  (29)  bears  a  very  similar  form  as  its  coun¬ 
terpart  for  the  optimum  energy  allocation  in  Eq.  (23).  In  fact, 
when  the  path  loss  exponent  v  =  2,  Eq.  (29)  is  essentially 
identical  to  Eq.  (23).  For  general  v  values,  however,  these 
two  relationships  are  quite  different.  Such  a  discrepancy  is 
actually  very  reasonable,  because  Eqs.  (23)  and  (29)  result 
from  two  distinct  optimization  problems:  the  former  is  ob¬ 
tained  for  arbitrary  distances  Ds  r  and  l)r.d  under  a  total  en¬ 
ergy  constraint;  whereas  the  latter  is  obtained  for  prescribed 
£s  and  £r  under  a  total  distance  constraint.  With  the  SER 
bound  Pe  being  a  two-dimensional  function,  the  energy  and 


Fig.  6.  SER  comparison  between  relay  systems  with  and  without 
energy  optimization  (p  =  lOdB,  v  =  4,  DBPSK). 

location  optimizations  are  carried  out  on  uncorrelated  dimen¬ 
sions. 

For  general  L  values,  the  optimum  location  can  be  deter¬ 
mined  in  the  similar  manner  as  we  discussed  in  the  preceding 
subsection.  Essentially,  the  path  loss  exponent  v  renders  it 
impossible  to  derive  the  analytical  solution  to  the  optimum 
location  problem,  even  with  the  high  SNR  approximation. 
One  can  resort  to  the  numerical  search  using  the  SER  bound 
in  Proposition  1.  In  Fig.  5,  the  optimum  distances  obtained 
from  the  numerical  search  and  the  simulations  are  compared 
for  different  v  values,  at  total  SNR  p  =  lOdB  and  with  L  1 
relay  node.  Notice  that,  as  its  counterpart  in  Fig.  4,  the  opti¬ 
mum  distance  allocation  is  linear  m£s/£  only  when  v  =  2. 

Summarizing,  we  established  the  optimum  energy  and  dis¬ 
tance  allocation  rules  which  minimize  the  SER  of  the  relay 
system.  The  analytical  solutions  were  derived  for  L  =  1,  and 
numerical  search  can  be  performed  for  L  >  2.  To  verify  the 
performance  improvement  resulted  from  the  optimization,  we 
will  present  extensive  comparisons,  generalizations  and  fur¬ 
ther  discussions  in  the  next  section. 

5.  SIMULATIONS  AND  DISCUSSIONS 

Next,  we  will  discuss  the  performance  of  the  DF  based  coop¬ 
erative  system  combined  with  differential  demodulation  and 
optimum  resource  allocation. 

To  verify  the  advantage  of  the  optimum  energy  allocation, 
we  plot  in  Fig.  6  the  SER  of  the  relay  system  with  and  without 
energy  optimization,  at  p  =  lOdB  and  with  various  numbers 
of  relays  L  =  (1, 2,  3, 4).  As  a  benchmark,  we  also  plot  the 
SER  of  a  direct  transmission  using  all  available  energy.  In  the 


Fig.  7.  SER  comparison  between  relay  systems  with  and  without 
location  optimization  (p  =  lOdB,  v  =  4,  DBPSK). 

system  without  energy  optimization,  a  uniform  energy  alloca¬ 
tion  is  employed;  that  is,  ps  =  pr  =  p/(L  +  1)  at  any  Ds  r. 
From  Fig.  6,  we  observe  that:  i)  the  un-optimized  system 
may  even  under-perform  the  direct  transmission,  regardless 
of  the  number  of  relays  L;  ii)  as  L  increases,  the  SER  perfor¬ 
mance  can  get  even  worse  unless  the  energy  optimization  is 
performed;  iii)  the  energy-optimized  system  universally  out¬ 
performs  the  direct  transmission  with  the  same  total  energy; 
and  iv)  the  energy-optimized  system  universally  outperforms 
the  un-optimized  system.  These  observations  confirm  our  dis¬ 
cussions  in  preceding  sections. 

Interestingly,  notice  that  the  minima  of  the  energy-optimized 
SER  curves  almost  coincide  with  those  of  the  un-optimized 
ones.  This  implies  that  the  near-optimum  SER  can  be  achieved 
even  with  the  uniform  energy  allocation  across  the  source  and 
relay  nodes,  provided  that  the  relay  location  is  carefully  se¬ 
lected. 

In  Fig.  7,  we  verify  the  advantage  of  the  optimum  dis¬ 
tance  allocation  by  comparing  the  SER  with  and  without  lo¬ 
cation  optimization.  Again,  the  SER  of  a  direct  transmis¬ 
sion  using  all  available  energy  is  plotted  as  the  benchmark. 
In  the  system  without  location  optimization,  the  relays  are 
placed  at  the  midpoint  of  the  source-destination  link,  as  sug¬ 
gested  in  [9,  13].  Similar  to  the  energy  optimization  case. 
Fig.  7  confirms  the  advantages  of  the  location  optimization: 
i)  the  un-optimized  system  may  even  under-perform  the  di¬ 
rect  transmission,  regardless  of  the  number  of  relays  L\  ii)  the 
location-optimized  system  universally  outperforms  the  direct 
transmission  with  the  same  total  energy;  and  iii)  the  location- 
optimized  system  universally  outperforms  the  un-optimized 
system.  Different  from  the  energy  optimization  case,  how¬ 
ever,  as  L  increases,  the  SER  performance  always  improves 
even  without  any  location  optimization. 


The  curves  in  Fig.  7  also  exhibit  more  flatness  compared 
with  the  ones  in  Fig.  6.  This  implies  that  the  system  SER 
is  more  sensitive  to  the  location  distribution  than  to  the  en¬ 
ergy  distribution.  In  addition,  the  minima  of  the  location- 
optimized  SER  curves  are  far  from  those  of  the  un-optimized 
ones,  except  for  the  L  =  1  case  (see  Fig.  7).  This  indicates 
that  placing  the  relay  nodes  at  the  midpoint  cannot  achieve 
the  minimum  SER  even  with  careful  allocation  of  the  source 
and  relay  energies,  for  any  L  >  1.  Notice  that,  even  with  opti¬ 
mum  energy  allocations  in  Fig.  7,  the  midpoint  relay  location 
still  cannot  achieve  the  minimum  SER,  unless  L  =  1.  This  is 
to  be  distinguished  from  the  uniform  energy  case  depicted  in 
Fig.  6,  as  well  as  from  the  coherent  relay  systems  in  [9]. 

6.  CONCLUSIONS 

In  this  paper,  we  investigated  the  optimum  energy  dis¬ 
tribution  and  optimum  location  of  relays  in  a  system  with 
arbitrary  number  of  relays  employing  differential  demodula¬ 
tion.  The  base  of  our  optimization  is  an  upper  bound  on  the 
average  SER,  which  we  derived  for  the  decode-and-forward 
cooperative  protocol.  Our  simulations  and  numerical  exam¬ 
ples  confirm  that  both  the  energy  and  location  optimizations 
provide  considerable  SER  advantages.  Without  optimization, 
the  system  with  more  relays  may  at  times  underperform  the 
system  with  less  relays  or  the  direct  transmission.  We  have 
also  shown  that  the  location  optimization  may  be  more  crit¬ 
ical  than  the  energy  optimization.  In  other  words,  the  dif¬ 
ferential  relay  system  with  uniform  energy  distribution  can 
achieve  near-optimum  SER  by  appropriately  choosing  the  re¬ 
lay  location;  while  a  system  with  relays  sitting  at  the  midpoint 
between  the  source  and  the  destination  cannot  approach  the 
optimum  SER  even  with  optimized  energy  distribution. 
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